Gamma-Ray Bursts Black hole accretion disks as a site for the z/p-process 
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We study proton rich nucleosynthesis in windlike outflows from gamma-ray bursts accretion disks 
with the aim to determine if such outflows are a site of the fp-process. The efficacy of this i/p-process 
depends on thermodynamic and hydrodynamic factors. We discuss the importance of the entropy of 
the material, the outflow rate, the initial ejection point and accretion rate of the disk. In some cases 
the vp-pTocess pushes the nucleosynthesis out to ^ ~ 100 and produces light p-nuclei. However, 
even when these nuclei are not produced, neutrino induced interactions can signiflcantly alter the 
abundance pattern and cannot be neglected. 

PACS numbers: 26.30.-k, 25.30.Pt, 97.60.Bw 



I. MOTIVATION AND METHODS 

Burbidge, Burbidge, Hoyle and Fowler [l[ and indepen- 
dently Cameron Q have proposed three major nucleosyn- 
thesis processes - the r-, s-, and p-process - to make nuclei 
heavier than iron in stars. Since this pioneering work im- 
pressive progress has been achieved in the understanding 
of these processes [3| . Nevertheless important questions 
still remain unanswered like the definite astrophysical site 
of the r- process [4] and the mechanism responsible for the 
production of light p-nuclei 5] . 

Neutrino-driven winds from core-collapse supernovae 
are considered as a site for the production of light p- 
nuclei. In this scenario, light p-nuclei can be produced 
by different mechanisms. Slightly neutron rich outflows 
with electron fractions Ye ~ 0.48 result in the production 
of Hght p-nuclei with A < 92, i.e. ^^Se, ^^Kr ^^Sr, and 
^^Mo The production of some of these nuclei and 
in particular ^^Mo is enhanced by neutrino captures on 
nuclei with N fa 50 [3]. Light-p nuclei, including ^"^Mo 
and ^^'^®Ru, can be produced in proton-rich outflows via 
the Z/'p-process [8l-[l0l|. It is this latter mechanism that 
we focus on in this paper. It operates in proton-rich en- 
vironments with high neutrino and antineutrino fluxes 
as are found in the early ejected matter in core-collapse 
supernovae. Here the matter is ejected from the surface 
of the nascent neutron star as free nucleons. The com- 
petition of neutrino captures on neutrons and antineu- 
trino captures on protons drives the matter proton-rich 
as both neutrino types have rather similar luminosities 
and the average antineutrino energy is not large enough 
compared to the neutrino energy to compensate for the 
difference in reaction Q-values. Upon reaching cooler re- 
gions, i.e. with increasing distance from the neutron star 
surface, the nucleons assemble in nuclei and, without fur- 
ther neutrino reactions, the proton-rich matter freezes 
out with a significant production of N = Z nuclei like 
^^Ni and ^''Ge and some free protons left. However, an- 
tineutrino captures on these protons ensure a significant 



presence of free neutrons which can be captured on the 
N — Z nuclei via {n,p) and {n,j) reactions allowing for 
matter flow beyond ^^Ni and ^^Ge which otherwise with 
their long halflives against proton capture and beta decay 
could not be overcome during the dynamical timescale 
of supernova nucleosynthesis (a few seconds). Different 
studies of the J/p-process |8l-[lO| have shown that it can 
produce light p-nuclides in abundances which might be 
sufficient to explain their observed solar values. It is ex- 
pected that the ejection of some proton-rich matter in the 
presence of strong antineutrino fluxes is a general feature 
of core-collapse supernovae making the occurrence of the 
i/p-process a general supernova phenomenon. 

Recent studies indicate that quite similar physical con- 
ditions, like those which give rise to the t/p-process in su- 
pernovae, also occur in the windlike outffows from the ac- 
cretion disk surrounding a black hole featuring a gamma- 
ray burst. Considerable observational evidence now con- 
nects long duration gamma ray bursts with the collapse 
of massive stars. McFadyen and Woosley [Tl| pioneered 
a currently widely used model for these events: a col- 
lapse of massive star with M > 25 Mq (Mq denoting 
our sun's mass) fails to produce a standard core collapse 
supernovae with a neutron star at the center, and in- 
stead produces an accretion disk surrounding a black 
hole. Here again, the temperature in the disk is suffi- 
ciently high so that matter is ejected as free nucleons. 
Furthermore, the ejection occurs in the presence of quite 
sizable neutrino fluxes Depending on the matter ac- 
cretion rate of the black hole disk, neutrinos and antineu- 
trinos can get trapped defining neutrino surfaces from 
which they decouple from the disk matter [l^, [l3 ■ It is 
found that disks with moderate accretion rates (around 
1 Mq s^^) have neutrino and antineutrino fluxes which 
drive the ejected matter proton-rich making the occur- 
rence of the lyp-pmcess an exciting possibility [l3| . Lack- 
ing detailed (magneto)hydrodynamic models of the ejec- 
tion, the amount of proton-rich ejecta can be estimated to 
reach up to 0.01 Mq according to the analytical outflow 
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model of ref. [l5| and typical disk parameters. A broad 
survey of nucleosynthesis from GRB accretion disks indi- 
cated that some p-process elements may be formed this 
way [165 , but a study of the j/p-process mechanism has 
not been undertaken. It is the aim of this paper to study 
this intriguing possibility. Our study is based on the 
temperature, density, electron-to-nucleon ratios profiles 
and dynamical expansion times of the windlike outflows 
as described by McLaughlin and Surman [l^ [l^ . The 
respective wind trajectories are then coupled to an ex- 
tensive nuclear network. 



II. MODEL 

The disks studied in this work are all based 
on the model of Neutrino-Dominated-Accretion-Flows, 
"NDAF's" Two NDAF-based disk models are used 
and results compared. The models labeled A1-A6 are 
from a disk model by T. DiMatteo, R. Perna and R. 
Narayan [l^ ; this model was the first to incorporate the 
effects of neutrino trapping. The models labeled B1-B5 
are from a fully relativistic disk model by W.-X. Chen 
and A.M. Beloborodov [I^ that additionally incorpo- 
rates improved microphysics, including the influence of 
electron degeneracy and an evolving neutron-to-proton 
ratio. For a detailed description of neutrino-cooled ac- 
cretion disks, see [l^, [13 ■ Radial proflles for the 
disk density p, temperature T, disk scale height H, nu- 
clear composition and neutrino and antineutrino lumi- 
nosity Li/^p are found to depend strongly on the accretion 
disk parameters, i.e., the disk accretion rate M, viscosity 
parameter a, black hole spin a and accretion disk mass 
M . The matter in the disk has to decrease its internal 
energy to be accreted into the black hole. This cooling of 
the disk proceeds either through advection or neutrino 
emission [19]. Cooling through neutrino emission only 
becomes important in regions where the temperature and 
density are sufficiently high for the below reactions to oc- 
cur with significant rates: 

+ ^ i^e + (la) 
e~-)-p ^ Ve + n (lb) 
+ n ^ Ve+p (Ic) 

These neutrinos emitted from the disk interact with 
material ejected from the disk, infiuencing the subse- 
quent nucleosynthesis in the outfiow. The neutrino and 
antineutrino fiuxes above the disk are calculated as in 
Surman and McLaughlin [l3| . Where the disk is optically 
thin to neutrinos and antineutrinos, the emitted neutrino 
fiuxes are determined from the rates of the above reac- 
tions. In the inner regions of the disk, first the neutrinos 
and then the antineutrinos become trapped. Here we 
determine the surfaces at which electron neutrinos and 
antineutrinos decouple vertically from the accretion disk 
and then use the local disk temperatures at these places 




FIG. 1: (Color online) Scheme of the geometry in core collapse 
supernovae and GRB accretion disks. "PNS" is the nascent 
proto-neutron star. Neutrinos and matter in the early ejecta 
are emitted radially. "BH" denotes the black hole, the central 
engine in a GRB accretion disk. The outflows leave the disk 
vertically. 

to adjust the temperatures of the neutrino and antineu- 
trino fluxes leaving the disk. 

The evolutionary path of mass elements in an outflow 
is given in a trajectory. The initial thermodynamic para- 
meters of the matter in the outflow are set by the char- 
acteristics of the disk, i.e. vary according to the disk 
model used. Matter in the outflows leaves the disk in 
perpendicular direction to the disk. Close to the disk a 
treatment in cylindrical symmetry is convenient. With 
increasing distance to the disk a transition from cylindri- 
cal to spherical symmetry is appropriate. 

The present calculations are based on the simulations 
of the windlike outflows of Surman and McLaughlin 
[T3 |. Using the method outlined in [ij] we construct 
parametrized trajectories which describe the ejection of 
mass elements from the disk. Due to the rather high 
temperatures the matter is originally ejected as free neu- 
trons and protons and the proton-to-neutron ratio is set 
by weak interactions due to the strong neutrino and an- 
tineutrino fluxes which occur rather close to the neutrino 
surfaces. Upon reaching cooler regions at distances no- 
ticeably further out than these surfaces nucleons can re- 
assemble to nuclei. This nucleosynthesis process in the 
matter outflow we describe by an extensive nucleosyn- 
thesis network subject to the thermodynamical condi- 
tions of the outflowing matter. 

We assume the outflows to be adiabatic deflned by the 
entropy per baryon S in units of k^. A second crucial 
parameter is the velocity of the ejected matter as a func- 
tion of distance R from the black hole in the center of 
the disk. This velocity u is parameterized as 

l"ht'oo(l-^y. (2) 

where i?o is the radius at which the outflow leaves the 
accretion disk (see Figure [1]), Voo is the asymptotic ve- 
locity of the matter (assumed to be 3 x 10^ km/s) and 
the parameter /3 determines the acceleration of the out- 
flow and therefore is referred to as "acceleration param- 
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eter" [T3|- The parameter j3 is iraportant as it defines 
the time the outflow is subjected to weak interactions in 
the strong neutrino fluxes. Large values of j3 correspond 
to slowly accelerating outflows with rather long neutrino 
interaction times, while small /3 values deflne outflows 
with fast accelerations and hence less time for neutrino 
interactions. 

As pointed out in [3] the electron fraction Yg of the 
ejected matter depends also critically on the rate by 
which the black hole accretes matter to the disk. Typical 
values are in the range 0.01-10 s~^ [TJ]. Disks with 
smaller accretion rates of order 0.1 s~^ and 1 M0 
s~^, as explored in our current work, are found to have 
only small regions of trapped neutrinos in the disk and 
even smaller ones for antineutrinos. Under such condi- 
tions there are relatively few antineutrinos interacting 
with the matter outflow, while neutrino captures on neu- 
trons are still sizable. For such low accretion rates the 
weak interaction will drive the outflowing matter proton- 
rich [l3| with free protons available making it a tempting 
site for j/p-process nucleosynthesis. It is exactly these 
conditions which we will study in the following. 

We note that for even smaller accretion rates neither 
neutrino nor antineutrino trapping occurs and the free 
proton fraction in the ejected matter is not significant. 
Such an environment is not favorable for j/p-process nu- 
cleosynthesis. 

Our nucleosynthesis network considers 3347 nuclei cov- 
ering the nuclear chart from protons and neutrons to 
^^^Eu. The appropriate reaction rates among these nu- 
clei mediated by the strong and electroma gne tic inter- 
action are taken from the reaction library '20']. Weak- 
interaction rates for nuclei are adopted from Zinner [2ll | . 
However, we find that only neutrino reactions with free 
protons and neutrons have influence on the nucleosyn- 
thesis process. The network has been numerically solved 
as outlined in fHlli. 



III. RESULTS 

As outlined above we expect that the nucleosynthesis 
in the accretion disk outflows depends on the parameter 
values for the mass accretion rate M, the acceleration pa- 
rameter /3, the radius i?o at which the matter decouples 
from the disk and the entropy in the outflow S. Further- 
more, the proton-to-neutron ratio of the ejected matter 
is crucial for the nucleosynthesis. As this ratio is set by 
weak interaction processes close to the neutrino surfaces 
it is expected to be different for the neutrino-dominated 
accretion disk models of DiMatteo et al. [18] and of Chen 
and Beloborodov [l^ which predict quite noticeable dif- 
ferences in the relative neutrino and anti-neutrino fluxes 
and spectra. Hence we have performed a set of nucle- 
osythesis calculations for both models varying the para- 
meters M, /3, i?o and S. The chosen parameter values 
for the various models which we studied in details are 
defined in Table HI where the labels 'A' and 'B' refer to 
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(km) 




(kB) 


Al & Bl 


1 


100 


2.5 


30 


A2 & B2 


1 


100 


0.8 


30 


A3 & B3 


1 


100 


2.5 


15 


A4 & B4 


1 


100 


2.5 


50 


A5 


1 


250 


2.5 


50 


B5 


1 


50 


2.5 


50 


A6 


0.1 


100 


2.5 


30 



TABLE I: Parameter values for the various studied models. 
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FIG. 2: (Color online) Model Al: Evolution of the cap- 
ture rates for positrons (solid line) and neutrinos (dashed) 
on neutrons and for electrons (dashed-dotted) and antineu- 
trinos (dashed-dotted-dotted) on protons as a function of the 
matter temperature. Notice that the rate for positron cap- 
ture on neutrons includes also the neutron beta decay rate. 
The evolution of the electron fraction, Ye is also shown (right 
y-scale). Temperature is measured from the moment matter 
decouples from the disk where it is still strongly neutron-rich. 



the disk models of [18j and [19|, respectively. 

We start with nucleosynthesis studies for the outflows 
from the disk model of DiMatteo et al. [li|. The Model 
Al with the modest accretion rate M = 1 Mq and 
relatively slow outflow velocities corresponds to condi- 
tions for which one expects the outflows to be proton- 
rich as at the radius Rq at which matter is expelled from 
the disk (100 km) the disk is dense and hot enough for 
neutrino producing reactions (Eq. [T]) to occur with signi- 
flcant rates. Figure [2l shows the time evolution of the 
electron fraction Yg and the capture rates of positrons 
and neutrinos on neutrons and electrons and antineu- 
trinos on protons. Time is measured from the moment 
matter decouples from the disk where it is still strongly 
neutron-rich and consists of free protons and neutrons. 
Their ratio, however, is changed due to fast electron- 
positron captures and interactions with the large neu- 
trino and antineutrino fluxes. Importantly the neutrino 
fluxes are about one order of magnitude larger than those 
for antineutrinos, while the average energy of antineutri- 
nos is around 2 MeV larger than the one of neutrinos. As 
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a consequence neutrino captures on neutrons dominate 
and drive the composition proton-rich and a peak value 
of Ye ^ 0.8 is reached. The increase in Yg is stopped, 
once alpha particles form which use up all available neu- 
trons. Due to the high particle thresholds of '*He neutrons 
are then protected against neutrino interactions, while 
the remaining free protons are still subject to strong an- 
tineutrino fluxes. Continuing antineutrino captures on 
protons supply more neutrons, which, in an "inverse" 
alpha effect [3, [13, [1^ are combined with additional pro- 
tons into more *He. As a consequence of these antineu- 
trino captures the Ye value decreases. At T ~ 3 GK 
the weak interaction rates have decreased sufficiently as 
the neutrino fluxes decrease with distance and the mat- 
ter composition freezes out with a Ye value around 0.72. 
Starting with the triple-alpha reaction a network of re- 
actions (called alpha process [1^ [13]) produces heavier 
nuclei, which can become the seed for additional nucleo- 
synthesis, and reduces the abundance of *He. For condi- 
tions with Ye > 0.5, the seeds consist mainly oi N = Z 
nuclei with multi-alpha structures like ^^Ni, ^°Zn, ^'^Ge, 
with an abundance of free protons Yp — 2Ye — f 3 . The 
long halflives of these seed nuclei against beta decay and 
proton capture would stop the nucleosynthesis flow, if it 
were not for the continuous supply of free neutrons pro- 
duced by antineutrino captures on the protons. This is 
the key element of the t'p-process. Indeed a sequence of 
{n,p) and (p, 7) reactions allows for synthesis of elements 
upto the mass range ^ ~ 80 — 100 for the conditions of 
model Al, as can be seen in Fig. [3l As a comparison 
we also show the abundance distribution obtained with 
the same nucleosynthesis network, however, switching off 
the neutrino and, importantly, antineutrino capture reac- 
tions once the freeze-out value of Ye is reached. Strikingly 
the matter flow stops at the alpha seed nuclei with long 
halflives (e.g. ^^Ni, ^°Zn..) as in this scenario no free 
neutrons are available to carry the matter flow to larger 
mass numbers as it is guaranteed in the j/p-process due 
to the late-time neutrons produced by antineutrino cap- 
tures on protons. 

We have seen that model Al allows for a substantial 
fp-process to occur. The outflow model A2 is identi- 
cal to model Al, except that it has a faster acceleration 
(smaller /? value). Thus, although the neutrino and an- 
tineutrino luminosities are the same as in model Al, the 
ejected matter has shorter time to interact with neutrinos 
and consequently the neutrino fluence is much smaller. 
While in model Al the neutrino fluence in the temper- 
ature range 3-1 GK is 4 x 10^^ cm~^ (a value similar 
to the one used in up-process nucleosynthesis studies in 
supernova environments [1, [§]) the value for model A2 
is only 10^^ cm~^. Consequently, model A2 is expected 
to show only a rather weak i/p-process which is indeed 
borne out by our nucleosynthesis calculation, as can be 
seen in Fig. |4l 

Another important parameter, which affects the mat- 
ter composition and the subsequent nucleosynthesis, is 
the entropy S. In general the larger the entropy, the 




FIG. 3: Final isotopic abundances (solid circles) for the nu- 
cleosynthesis calculations in disk model Al relative to solar 
abundances [28|. The empty circles show the abundances 
obtained with the same nucleosynthesis network, however, 
switching off neutrino and anti-neutrino capture reactions 
once the freeze-out value for Ye is reached. 
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FIG. 4: Same as Fig. H but for model A2. 
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FIG. 5: Same as Fig. [1 but for model A3. 
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FIG. 6: Same as Fig. El but for model A4. 
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FIG. 7: Same as Fig. El but for model A5. 



larger the fraction of free protons which are available per 
seed nucleus. To study the influence of entropy on the 
final nucleosynthesis abundances we have repeated the 
parameter study of model Al, however, replacing the en- 
tropy of the ejected matter by a smaller value (15 fcs, 
model A3) and by a larger value (50 ks, model A4) 
than used in model Al (30 fcs). We mention again that 
the neutrino fluxes are the same in models Al, A3, and 
A4. The fact that for model A3 (15 ks), with lower 
entropy, there are less free protons available per seed 
nuclei translates into a smaller supply of neutrons pro- 
duced by antineutrino captures and a less pronounced 
Z/'p-process as observed for model Al. The behavior is 
opposite for model A4 (50 fc^), where the larger entropy 
reults in more free protons per seed nucleus, and con- 
sequently more free neutrons per seed, than model Al. 
These observations are confirmed by the abundance dis- 
tributions shown in Figs. [5] and [6| Model A4 (50 ks) 
indeed shows a substantial amount of elements in the 
mass range A ~ 60 — 100, while model A3 (15 ks) pro- 
duces only nuclides up to mass A w 80 in a noticeably 
weaker z/p-process. 

As mentioned before, the disks have regions of differ- 



n(e+,v Jp + n(P")p 
n(Ve,e")p 
p(e",Ve)n 
p(Ve,e+)n 




0.65 



0.60 



0.55 



0.50 



0.45 



Temperatue (GK) 
FIG. 8: (Color online) Same as Fig. [H but for model A6. 
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FIG. 9: Same as Fig. El but for model A6. 



ent nuclear composition, density and temperature. The 
weak interaction rates in the outflows depend on the neu- 
trino and antineutrino fluxes originating from the disk. 
The flux magnitudes in turn reflect disk conditions at 
certain radii. Thus the radius i?0: at which matter is 
released from the disk, affects t'p-process nucleosynthesis 
as weak interaction rates in the outflow depend on the 
magnitude of the neutrino and antineutrino fluxes. In 
disks with comparable accretion rate at larger radii neu- 
trinos will drive the composition less proton-rich as due 
to disk composition the neutrino fluxes are smaller. At 
smaller radii neutrinos will drive the composition more 
proton-rich and neutrino fluxes are larger. 

To demonstrate this quantitatively we have repeated 
the nucleosynthesis calculation of model A4, however, re- 
placing the radius at which matter is released to a larger 
value of Rq — 250 km (model A5). As expected, a weaker 
fp-process occurs in model A5 (Fig. [7]). 

Model A6 is the same as model Al, except that the 
matter accretion rate is significantly lower, M = 0.1 M0 
s~^ rather than 1 Mp^ as assumed in all other models. 
As is explained in [ij this change in accretion rate has 
significant effects on the (anti-)neutrino surfaces. Conse- 
quently it also strongly affects the neutrino and antineu- 
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FIG. 10: (Color online) Same as Fig.[2l but for model Bl. 



trino fluxes which now are quite similar rather than domi- 
nated by a much stronger neutrino than antineutrino flux 
as encountered in the other models. The fluxes are in 
overall a factor of magnitude smaller than in models with 
higher accretion rate. Although the average antineutrino 
energy is slightly larger than the average neutrino en- 
ergy, it is not the neutrino capture reactions which drive 
the composition slightly proton-rich in this model (see 
Fig. [5]) but positron captures on free neutrons. The 
rates for the discussed reactions (Eq. [1]) depend on the 
temperature and density of the disk. In fact, in model 
A6 a sufficient amount of positrons is produced which 
drives reaction ([Tc]) to the right and makes the ejected 
matter proton-rich, reaching Ye values up to 0.55. As 
is explained for model Al, once nuclei can be formed 
all neutrons are blocked in "^He with some free protons 
left. The formation of '^He occurs before weak freeze 
out. However, this value for Yf, is significantly smaller 
than in model Al (Ye ~ 0.72). This fact together with 
the smaller antineutrino fluxes results in a smaller pro- 
duction of neutrons per seed nuclei and a rather weak 
i/p-process, as is demonstrated in Fig. [S] 

We now switch to the nucleosynthesis studies for 
the outflows from the disk model of Chen and Be- 
loborodov [l^. This model predicts neutrino fluxes and 
spectra which are quite distinct from those of the disk 
model of DiMatteo et al. [l^ • The differing treatments 
of relativity and the microphysics of the two disk mod- 
els results in markedly different disk temperature and 
density profiles, particularly in the innermost regions 
with the greatest neutrino emission. In the DiMatteo 
et al. disk model, the temperature and density rise 
more steeply with decreasing radius and the innermost 
regions are significantly hotter and denser than in the 
Chen and Beloborodov disk model, resulting in corre- 
spondingly higher neutrino emission. 

To explore the impact of these differences on the nu- 
cleosynthesis we have performed nucleosynthesis studies 
for the Chen and Beloborodov model (called models Bl- 
B4) using the same parameters for the mass accretion 
rate, the acceleration parameter, the decoupling radius. 




Mass Number A 
FIG. 11: Same as Fig.E] but for model Bl. 



and the entropy as in models Al- A4 above (see Table HI . 
Fig. [ini is the equivalent to Fig. [2] clearly showing im- 
pact of the different neutrino fluxes on the proton-to- 
neutron ratio of the ejected matter. In contrast to model 
Al, in which the neutrino luminosity exceeds the one for 
anti-neutrinos, by about an order of magnitude, the neu- 
trino and anti-neutrino fluxes are rather similar in model 
Bl. As a consequence also the neutrino and anti- neutrino 
capture rates on neutrons and protons, respectively, differ 
signiflcantly less than in model Al. Moreover, neutrino 
and anti-neutrino capture rates are both smaller than the 
rates for positron and electron captures during the flrst 
second after matter has decoupled from the disk. This 
situation is similar to model A6, and again it is the slight 
dominance of positron over electron captures which drive 
the matter proton-rich at flrst. There is only a little 
time period in model Bl for which neutrino and anti- 
neutrino captures dominate the inverse weak-interaction 
processes. As the anti-neutrino capture rate is slightly 
larger by a factor of order 2, the ejected matter is fur- 
ther driven proton-rich during this short period. The 
flnal proton-to-neutron ratio Yg = 0.56 is much smaller 
than in model Al (Ye = 0.72). We also observe an in- 
verse alpha-effect in model Bl, but it is milder than in 
model Al. Since Ye > 0.5, the z/p-process can also oper- 
ate in model Bl. However, its effectiveness is noticeably 
less than observed in model Al, as less free protons are 
available with a reduced antineutrino flux (see Fig. [TT]) . 
While nuclides beyond A = 64 are being made, the i^p- 
process stops in model Bl already around A — 80 and 
produces also smaller abundances of nuclides in the mass 
range A = 60 — 80 than model Al. 

Model B2 is the same as model Bl, except that it has 
a smaller value for (3 (faster acceleration). Hence weak- 
interaction processes have less time to change the proton- 
to-neutron ratio. In fact, as can be seen in Fig. [T^ the 
flnal value for Ye is already been reached at a time of 0.1 
s after the matter decoupled from the disk. More impor- 
tantly for the nucleosynthesis, the value of Ye stays below 
0.5 in this model, i.e. the ejected matter is neutron-rich. 
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FIG. 12: (Color online) Same as Fig. [1 but for model B2. 
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FIG. 13: Same as FigO but for model B2. 



We also observe that positron and electron captures dom- 
inate over neutrino and anti-neutrino captures during the 
time period in which the Ye value of the ejected matter is 
set to its final value by weak interactions. With the mat- 
ter being slightly neutron-rich, obviously no z^p-process 
occurs. Fig. [T3l shows that the abundance distributions 
are nearly identical with and without consideration of 
neutrino and anti-neutrino reactions in the nuclear net- 
work after the final Ye value has been reached. These 
abundance distributions resemble those of an a-process 
for a given value of Fe < 0.5. The distinct differences in 
the abundance distributions of an alpha-rich freeze-out 
for proton-rich (Ye > 0^ and neutron-rich (Ye < 0.5) 
matter is explained in ^^.) We note that such an a- 
process operates also in the neutrino-driven wind scenario 
setting up the abundance distribution for the seed nuclei 
of an subsequent r-process. However, in our model B2 
the entropy is significantly smaller than required for suc- 
cessful r-process simulations in the neutrino-driven wind 
model. Due to this low entropy the neutron-to-seed ratio 
is too low in model B2 to allow for any r-process to occur. 
We note that the nuclear network considered in our nu- 
cleosynthesis studies is large enough to indicate potential 
r-process nucleosynthesis beyond the a-process. 
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FIG. 15: Same as Fig. El but for model B4. 



Models B3 and B4 are the same as model Bl, how- 
ever, assuming a smaller or larger value for the entropy 
(5' — ISfcs, model B3 and S = bOks, model B4, re- 
spectively). As explained above, we expect a more pro- 
nounced J^p-process to occur with growing entropy, as a 
larger amount of free protons (and hence neutrons af- 
ter anti-neutrino captures) is available once heavy nuclei 
are being formed in the ejected matter. In fact, this is 
born out by our nucleosynthesis studies for models B3 
and B4 (see Figs. rHl and [TSl) . if compared to model Bl. 
It is also interesting to compare the abundance results 
of models B3 and B4 to the equivalent studies of the 
DiMatteo et at disk model (models A3 and A4). In 
both cases, the stronger i^p-process nucleosynthesis oc- 
curs for the DiMatteo et al. disk models. The reason 
is the same as discussed in the comparison of models 
Al and Bl. The reduced enhacement of neutrino over 
anti-neutrino fluxes in the Chen and Beloborodov models 
translates into smaller Ye values than in the correspond- 
ing DiMatteo et al. models. The consequence are smaller 
proton-to-seed (or neutron-to-seed) ratios available once 
the i^p-process gets operable. 

Finally, in model B5 we explore the effect of reducing 
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FIG. 16: Same as Fig. El but for model B5. 



the radius Rq at which matter is released from the disk. 
This model is identical to B4, however, it uses a value for 
Rq of 50 km rather than 100 km. As discussed above, this 
change leads to an increase in the neutrino fluxes that will 
drive the composition more proton-rich and to a stronger 
^'p-process. This expectation is indeed borne out in our 
calculations (see fig. [TO)) , which shows noticeable strong 
production of elements with A > 75. In contrast to B4 
(see fig. [T5]) , model B5 synthesizes also elements in the 
mass range A = 90-100. 



IV. CONCLUSIONS 

Gamma-ray bursts are likely connected to an accretion 
disk surrounding a black hole. It has been found that 
the matter ejected from this disk can exhibit a rather 
rich range of nucleosynthesis fl6i\ dependent on the mat- 
ter accretion rate of the disk, the velocity and entropy 
of the outflow as well as the radius at which the mat- 
ter is released. The composition of the outflow can be 
either neutron-rich, giving range to r-process like nucleo- 
synthesis fl3| , or proton-rich. In this manuscript we have 
chosen the parameters characterizing the outflow such 
to focus on proton-rich nucleosynthesis and to investi- 
gate whether a up-piocess can occur under reasonable 
conditions in the windlike outflows from a GRB accre- 
tion disk. This is indeed confirmed by our extensive 
nuclear network simulations. We find that it is essen- 
tial to include neutrino capture interactions, particularly 
anti-neutrino capture on protons, when determining the 
nucleosynthetic outcomes for these environments. 

In our present calculations, including disk models 
with accretion rates from 0.1 Mq s~^ and 1 Mq s~^, 
we find that disk models with accretion rates of order 
1 Mq s~^ are particularly favorable for i/p-process nu- 
cleosynthesis as these models exhibit significantly larger 
neutrino fluxes than antineutrino fluxes driving the mat- 
ter proton-rich by a competition of neutrino captures on 
neutrons and antineutrino captures on protons. Our net- 



work calculations show that for reasonable parameters 
describing velocity, entropy and ejection radius of the 
outflow a ^'p-process can occur which produces substan- 
tial abundances of nuclides up to the A ^ 100 mass range. 
The effectiveness of the production in the mass range 
A r-^ 60-100 depends somewhat on the entropy of the 
outflow; where larger entropies support a larger amount 
of free protons per seed nucleus and hence give a stronger 
vp-pvocess as, by antineutrino capture, the free protons 
supply the neutrons for up-process nucleosynthesis. 

A crucial parameter for successful i^p-process nucleo- 
synthesis is the acceleration parameter /? of the ejected 
matter as it defines the time matter is subjected to the 
neutrino interactions. Obviously the bigger the accele- 
ration parameter, the more neutrino captures drive the 
matter proton-rich. Indeed by varying the acceleration 
parameter by a factor 3 we have demonstrated that the 
nucleosynthesis outcome can change from the occurrence 
of a strong i/p-process to no i^p-process nucleosynthesis. 
Slightly weaker (stronger) z/p-process nucleosynthesis is 
also observed if the matter is released at larger (smaller) 
radii from the disk due to reduced (enhanced) neutrino 
capture rates. 

It is also worth mentioning that our calculations show 
the occurrence of a process which in analogy to the al- 
pha effect in r-process nucleosynthesis 0, [24,[25| we like 
to call "inverse" alpha effect. If the formation of '^He, 
which in proton-rich environment locks up all neutrons, 
occurs already before weak interaction freeze out, then 
continuous antineutrino capture on protons supply more 
neutrons, which combine with protons to build more ^He. 
This inverse alpha effect reduces the proton-to-neutron 
ratio in proton-rich environments (while in neutron-rich 
conditions it increases Ye). The inverse alpha effect is 
particularly important for outfiows with rather slow ac- 
celeration parameters. 

Accretion disk models are currently under development 
by several groups, so we have compared nucleosynthesis 
for the outfiows of two different accretion-disk models. 
While we find the same basic conclusion: that the u-p 
process occurs, the details depend on the particular ac- 
cretion disk model used. The model of DiMatteo et al. 
predicts a noticeable excess of neutrino flux over anti- 
neutrino flux. This difference is much milder in the disk 
model of Chen and Beloborodov. As a consequence of 
these differences in relative fluxes, the ejected matter in 
the disk model of Ref. 1181 ha s always larger values 
than found for the model of [1^ . Obviously the occurence 
of an z/p-process is always more pronounced for the for- 
mer model if nucleosynthesis abundances are compared 
for the same values of mass accretion rate, matter accele- 
ration, decoupling radies, and entropy. With the models 
A2 and B2 we have presented examples, where one disk 
model predicts the occurence of an z/p-process, while, 
for the same parameter values, for the other [l^ matter 
is ejected with a Y^ value less than 0.5 translating into 
an abundance distribution known from a-rich freeze-out 
for neutron-rich matter. 
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In summary, within our parameter studies we have 
shown that outflow from accretion disks surrounding a 
black hole can be the site of t'p-process nucleosynthesis. 
More quantitative studies of this exciting perspective 
have to wait until complete hydrodynamical simulations 
for the outflows from the accretion disks become avail- 
able. 
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